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We study the uniform magnetic susceptibility and spin- 
gap for recently proposed Heisenberg model Hamiltonians for 
CaV40g based on the orbital ordering scenario of Marini and 
Khomskii and the LDA calculations of Pickett. We argue that 
the experimentally observed uniform susceptibility data is in- 
consistent with the weakly coupled dimer picture of Marini 
and Khomskii. The model proposed by Pickett can, with ap- 
propriate choice of parameters, lead to an explanation for the 
the observed gap and uniform susceptibility. The resulting 
agreement with experiments is of a similar quality to previ- 
ously studied models. We argue that this new model is best 
distinguished from previous ones by neutron or Raman scat- 
tering experiments, via the location of the excitation mini- 
mum in the Brillouin zone and by the possible existence or 
non-existence of sharply defined singlet excitations. 



The observation of spin gaps in the calcium vana- 
dates, especially in the quasi-2D material CaV^g, has 
attracted much attention. A number of theoretical pa- 
pers have explored the compatibility of the experimen- 
tal data with various model spin Hamiltonians. Ignoring 
some out-of-plane distortions, a given layer of these ma- 
terials forms a one-fifth depleted square lattice with a 
four-spin unit cell (see Fig. |l|). It is by now well estab- 
lished that a uniformly nearest-neighbor coupled, one- 
fifth depleted square-lattice Heisenberg model has long 
range magnetic order Ha Thus the spin gap has been ar- 
gued to arise either (i) from non-uniform couplings which 
can lead to weakly coupled plaquettes or (ii) from frus- 
trating second-neighbor interactions, which destabilize 
the Neel order. Either situation leads to a "plaquette 
RVB" (PRVB) ground state, and lies in the same T = 
phase as a Hamiltonian which contains only intraplaque- 
tte nearest-neighbor interactions. Based on extensive 
comparison of the experimental susceptibility and spin 
gap with theoretical models, it was argued by Gelfand 
et alau that the experimental data is best supported by 
the latter scenario, with the second-neighbor interactions 
roughly half the nearest-neighbor interactions. 

In a recent paper, Marini and KhomskiiQ suggest that 
the magnetic properties of the vanadates CaV„02n+i 
(with n = 2, 3, 4) are determined principally by orbital 
ordering and associated structural distortions. These 
considerations lead to rather different spin Hamiltonk 
ans than those considered by most authors previouslypel 
who neglected the orbital degrees of freedom. To leading 



order Marini and Khomskii find the system to be a set 
of weakly coupled dimers. _ 

In another recent paper, PickettEl has used LDA calcu- 
lations to study the magnetic properties of CaV^g. His 
most intersting conclusion is that the relevant orbitals are 
d x 2_ y 2 — different from the ones considered by Marini 
and Khomskii. Furthermore, he argues that out-of-plane 
distortions in the arrangement of the vanadium atoms 
are quite important. The one-fifth depleted square lattice 
of vanadium atoms is actually made up of two (planar) 
layers, slightly above and below one a6-plane, and each 
forming a square lattice of "meta-plaquettes." These 
meta-plaquettes are to be distinguished from the plaque- 
ttes considered before by Ueda et alu and other authors 
in the "plaquette RVB" (PRVB) scenario. The plaque- 
ttes involve nearest-neighbor spins, two of them lying in 
the upper layer and two in the lower layer. In this case, a 
significant second-neighbor interaction is expected within 
the plaquettes. In contrast, the meta-plaquettes of Pick- 
ett lie cither entirely in the upper or entirely in the lower 
layer. Pickett argues that although the spins within a 
meta-plaquette are further apart than spins between two 
neighboring meta-plaquettes, various quantum-chemical 
arguments conspire to make the interactions within meta- 
plaquettes the strongest. He finds that to leading ap- 
proximation the magnetic system should be regarded as 
decoupled meta-plaquettes. One important point to note 
is that within a meta-plaquette the interactions are ex- 
pected to be primarily between neighboring spins only. 

The aim of this paper is to examine the consistency 
of the above two scenarios with existing and prospec- 
tive experimental data on CaVjOg. Neither Marini and 
Khomskii nor Pickett provide quantitative estimates for 
all the relevant exchange constants. Within their scenar- 
ios one might easily consider a four parameter Hamilto- 
nian in each case. In the interest of simplicity and be- 
cause our main conclusions are unaffected by the number 
of parameters, we have considered only a two-parameter 
Hamiltonian in both cases. 

First we consider the proposal of Marini and Khomskii. 
They suggest that the spins in CaV40g form antiferro- 
magnetic dimers with exchange interaction J approxi- 
mately equal to the observed triplet gap, a bit larger than 
100 K. If the inter-dimer interactions are neglected then 
the susceptibility follows \{T) oc T- 1 /(e J / T + 3). This 
form is shown as the dashed line in Fig. ^| and compared 
with the data of Taniguchi et alS3 (solid line). Here we 
have chosen J and the normalization for the susceptibil- 
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ity so that the different curves coincide at the maximum. 
The two curves are clearly inconsistent. Including the 
"corner" exchange constants (shown by thin solid lines 
in Fig. |l|), suggested to be the second most imporant by 
Marini and Khomskii, leads to an alternating spin-chain 
Hamiltonian 
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Here a = corresponds to the dimer limit. Accu- 
rate calculations fop-.this model are possible by finite- 
size diagonalization.Eil Once again we adjust the exchange 
constants and the g-factors to make the theoretical sus- 
ceptibility coincide with the experimental data at its 
peak. As shown in Fig. |^, increasing a to 0.5 still leaves 
one far from agreement with the data. It should be noted 
that much better fits are obtained within the PRVB 
scenariou. One can trace the inconsistency of the alter- 
nating chain model with the susceptibility data to the 
fact that the model's Curie constant is J(l + a)/4, which 
implies a much larger J value than what is needed to 
account for the gap. 

We now turn to the scenario presented by Pickett. 
Here the leading exchange constant J couples spins 
within a meta-plaquette. The second most significant 
interaction is that between nearest-neighbor vanadium 
spins. These are the nearest-neighbor plaquette and 
dimer interactions of Ueda et al. (shown by thick solid 
and dashed lines in Fig. [I]), and we take them to be equal 
to each other. They couple spins on neighboring meta- 
plaquettes with three terms in the Hamiltonian. Let 
these interactions be of strength AJ. The spin-gap for 
this model as a function of A is calculated via perturba- 
tion theory to be 

A/J= I - (2/3)A - 0.25347222A 2 - 0.048747749A 3 
-0.069181496A 4 - 0.030457074A 5 + 0(A 6 ) . (2) 

Dlog-Pade approximants suggest that this model orders 
magnetically at A c ~ 0.88, with a gap exponent approxi- 
mately 0.7, consistent with a d — 3 Heisenberg universal- 
ity class. The Curie- Weiss constant for the model is given 
by J(l+3A/2)/2. Thus the observed gap and Curie- Weiss 
constant can both be accounted for with J w 220K and 
A ~ 0.6. We note that these parameters will change 
if interactions between second neighbor metaplaquettes 
(that is, metaplaquettes in the same layer) are considered 
or if the nearest-neighbor dimer and plaquette interac- 
tions are unequal. However, the quality of the agreement 
with the high temperature susceptibility data is likely 
to be similar to that obtained within the PRVB scenario. 
The high temperature susceptibility is determined by the 
Curie- Weiss constant and the g-factor. Just as in previ- 
ous fits within the PRVB scenario, ^-values smaller than 
2 (roughly 1.7) will be needed to fit the data. 

The question remains as to whether one can exper- 
imentally distinguish a weakly coupled meta-plaquette 
scenario from a PRVB phase stabilized by frustration 



scenario. One approach is to measure the dispersion of 
triplet elementary excitations. In Pickett's model, the in- 
teractions between neighboring meta-plaquettes are two 
plaquette couplings and one dimer coupling. Assuming 
those two types of couplings to have equal strength, as 
in our derivation of Eq. (Q) , leads to a spectrum in which 
the gap minimum occurs at wavevector (0, 0). (Note that 
this corresponds to incipient two-sublattice order, not in- 
cipient ferromagnetism, due to the geometry of the meta- 
plaquettes.) The dimer couplings would need to be nearly 
twice as strong as the plaquette couplings in order for the 
gap minimum to occur elsewhere in the Brillouin zone. 
In the PRVB scenario considered earlier the gap mini- 
mum most likely lies at an incommensurate wavevector, 
or possibly at (7r,7r). Thus inelastic magnetic neutron 
scattering experiments may help to distinguish between 
these cases. Frustrating interactions within piaquettes 
can lead to sharply defined singlet excitationsjaiiia which 
are absent in Pickett's model. Thus Raman scattering 
can also help distinguish between the two scenarios. 
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FIG. 1. One- fifth depleted square lattice representing vana- 
dium spins associated with one afe-plane of CaV40g. The 
shaded circles represent vanadium atoms slightly above the 
plane and the open circles represent vanadium atoms slightly 
below the plane. The thick solid lines represent the elemen- 
tary plaquettes of Ueda et al. The thin dashed lines represent 
the meta-plaquettes of Pickett. In Marini and Khomskii's 
picture the largest interactions are the "dimer interactions" 
shown here by thick dashed lines and the second largest ones 
are the "corner" interactions, shown by thin solid lines. In 
Pickett's picture the largest interactions are the metaplaque- 
tte interactions shown by the thin dashed lines and the sec- 
ond most important interactions are the dimer and plaquette 
interactions shown by the thick dashed and solid lines respec- 
tively. 

FIG. 2. Uniform susceptibility for CaV409 compared with 
the alternating spin-chain model with a — (dimer limit) and 
a — 0.5. The exchange constant and the normalization for the 
susceptibility are scaled to make the curves coincide at their 
maxima. 
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